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ABSTRACT 

The Kepler mission was launched in 2009 and has discovered thousands of planet candidates. 
In a recent paper, [Esteves et aT] ( |2013| l found a periodic signal in the light curves of KOI-13 
and HAT-P-7, with a frequency triple the orbital frequency of a transiting planet. We found 
similar harmonics in many systems with a high occurrence rate. At this time, the origins of 
the signal are not entirely certain. 

We look carefully at the possibility of errors being introduced through our data process¬ 
ing routines but conclude that the signal is real. The harmonics on multiples of the orbital 
frequency are a result of non-sinusoidal periodic signals. We speculate on their origin and gen¬ 
erally caution that these harmonics could lead to wrong estimates of planet albedos, beaming 
mass estimates, and ellipsoidal variations. 
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1 INTRODUCTION 

To date, the Kepler mission has discovered over one thousand con¬ 
firmed planets and many more planet candidate^ The main goal of 
the mission was to discover Earth-like planets but many other sci¬ 
entific results have made use of this high quality dataset ( |Borucki| 
|2010| l. For example, |Barclay| ( |2012| > used the Kepler dataset to mea¬ 
sure ellipsoidal variations which can be used to derive masses and 
radii of both planets and stars. 

In a recent paper, [Esteves et al.| ( |2013| > found a mysterious 
third harmonic in the phase curves of Kepler-13b and HAT-P-7b. 
Until now, this signal had been detected in few planetary systems. 
Previously, tidal effects were the only known source of this effect 
( [Morris & Naftilan||1993^ . In our own investigation, we observe 
these higher order harmonics in a number of systems, including 
Kepler-13b, but tidal effects do not provide a sufficient explanation. 
The mostly likely source of these signals, is due to non-sinusoidal 
periodic variations in the light curves of these planets. 

The outline of the paper is as follows. First, we describe the 
methods and results in Sect.|^ In Sect.[3we discuss and dismiss the 
possibility of introducing this signal in our pipeline. The analysis 
of tidal effects are covered in Sect. [4] We summarize and conclude 
our results in Sect. [5] 


* caden.armstrong@mail.utoronto.ca 
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2 METHOD AND RESULTS 

We worked with publicly available data of Kepler light curves 
( [Borucki[2010^ . The planets selected were those with the confirmed 
tag in the Kepler Object of Interest (KOI) table from the Mikul- 
ski Archive for Space Telescopes (MASTQ We looked at the long 
cadence data set including quarters 0-17. The data was retrieved 
from the .fits files using the PyFits module, available in the AstroPy 
project ( [Astropy Collaboration et al.[2013^ . 

The Kepler team provides different data products for each 
system. Here, we investigate the PDCSAP_FLUX and SAP_FLUX 
columns. The SAP_FLUX column contains the light curve data pro¬ 
duced using simple aperture photometry which was produced by 
the Photometric Analysis module of the NASA Science Opera¬ 
tions Center pipeline. The column provides flux in electrons per 
second contained in the optimal aperture pixels, collected by the 
Kepler spacecraft. The PDCSAP_FLUX column contains the Photo¬ 
metric Analysis data after the Pre-Search Data Conditioning (PDC) 
algorithm has been applied. The PDC module applies a series of 
corrections to the data for known instrumental anomalies, along 
with unanticipated artifact^ Our analysis was run on both the 
PDCSAP_FLUX and SAP_FLUX columns with little to no effect on 
our final results. 

The flux for each planet was normalized to the average flux in 
each quarter. Normalizing to the average provided the simplest and 
most consistent method of combining the different quarters with¬ 
out introducing artificial signals. The values for planet period, tran- 

^ Updated March 1, 2013 

^ Kepler Data Characteristics Handbook. 31 May 2013. 
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Frequency (1/orbital period) Phase 


Figure 1. The left panel show the periodograms plotted in a log scale, with power in y-axis and the false alarm probabilities plotted as the dotted lines. The 
figures on the right are the folded light curves for illustrative purposes, showing the removed sections of data. The plot scaling is set to only show analyzed 
points, both the transits and outlying data points may extend outside the visible area of the plot. Blue data points were removed prior to creating the power 
spectra. The black line is the fitted three harmonics model, and the grey dots are binned data points. 


sit midpoint, and transit width were taken from the MAST Kepler 
catalogue. The separate light curves from each quarter were then 
merged together to create a single light curve. 

All data points falling within the primary and secondary transit 
of each planet were removed from the light curves. The width of the 
transits was increased by a factor of two to account for any transit 
timing variations, duration variations, or small errors in the cutting 
process. By comparison, the transit duration variation as measured 
in [Mazeh et al.| p013| > are only upwards of a few percent. Note 
that cutting out all transits with a longer transit duration reduces 
the signal to noise ratio but proved to be the simplest method of 
ensuring that the transits would not contribute to the signal (see 
also discussion in Sect. [^. The folded light curves are shown on 
the right panel of Fig. [^nd Fig. Data points in blue are those 
near the primary or secondary transit and have been cut. In addition 
to the transits, outlying data points were removed. 

To verify the period of the signal, we use a power spectrum. 
We turn to a Lomb-Scargle ( [Press & Rybicki|1989| > periodogram as 
the data points are unevenly spaced. Lomb-Scargle periodograms 


were created using the module available in the SciPy package 
( [Jones et al.|200l| l. The periodograms were sampled in an evenly 
spaced log scale from 10~® to 10® * orbital frequencies, with 2000 
sample points. The power spectra were plotted, and analyzed by 
eye, looking for any spikes at frequencies that correspond to mul¬ 
tiples of the orbital frequency. Only spectra that had a high peak 
compared to the overall noise were considered. This left us with 16 
systems with clearly identifiable signals with a high signal to noise 
ratio. Their properties are listed in Table The signal found by 
[Esteves et al.[ ( [2013) in HAT-P-7 was not recovered. 

The periodograms for these systems are plotted in the left pan¬ 
els of Fig.[T]and Fig.|^ They show the following general features: 


(i) A peak on the orbital frequency of the planet, attributed to 
Doppler boosting (also known as beaming, Rybicki & Lightman[ 
[1979[ l and planetary reflection ( [Faigler & Mazeh[2011[ |. 

(ii) A second peak on twice the orbital frequency, attributed to 
the ellipsoidal deformation of the star under the effects of the planet 
( [Barclay[2012[ l. 
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Frequency (1/orbital period) Phase 


Figure 2. Same as Fig.^ 


(iii) Additional harmonics at multiples of the orbital frequency 
of unknown origins. 

Of the 16 systems, the third harmonic is easily visible in 13 of them. 
In all of these systems, the periodograms showed narrow peaks 
exactly at the triple orbital frequency. In a few of these systems, 
similar peaks are found at higher orbital frequencies. The width of 
the peaks for these harmonics are of similar size compared to the 
known lower order harmonics. It is important to note that many of 
the systems showing the signal have been mentioned in the two pa¬ 
pers that look at transit timing variations of Kepler planets. The fol¬ 
lowing planets are mentioned in|Mazeh et al.|p013|l: KOI-823.01, 
KOI-1546.01, KOI-1459.01, and KOI-1285.01. The following are 
mentioned in |Szabo et allpOTS) : KOI-256.01, KOI-823.01. While 
not all of the transit timing variations were measurable in |Mazeh| 
|et al.|p013^ , they were still concluded to be significant. It is also 
interesting to note that in |Mazeh et al.|pOT3) , the third harmonic for 
KOT203.01 is visible in their Fig. 17 but not discussed any further. 
In our own analysis, the signal to noise ratio was not significant 
enough in KOI-203.01 for signals to be visible. 

In the MAST table for the Kepler Objects of Interest, the cho¬ 
sen planets were all of the CONFIRMED status, as of March 1st, 


2013. In |Muirhead et al.|{2013| l, the Kepler Object of Interest KOI- 
256 was concluded not to be a planetary system, but a mutually 
eclipsing binary consisting of a cool white dwarf and an active M3 
Dwarf. The expected transit is in fact the occultation of the white 
dwarf passing behind the M dwarf. The period for the white dwarf 
is unchanged from the value reported in the MAST table. In their 
analysis, higher order harmonics of the white dwarf are found in 
the light curve of the system. 


3 ARTIFACT OR REAL SIGNAL 

The first possibility of the signal’s origins resides in the misinter¬ 
pretation of the periodograms. It is easy to assume that the signal 
is produced as a combination of the known signals (beaming, el¬ 
lipsoidal variation, etc.). However, there is no way to combine the 
lower order periodic signals on the period of 1 orbital period and 

4 orbital period, to create the higher order harmonic. To check for 
these types of errors computationally, synthetic light curves were 
generated with the known signals, and analyzed with the same 
pipeline used to process the Kepler data. From these tests, no higher 
order harmonics could be recovered. 
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KOI Name Period [days] Depth [ppm] Planet radius [r®] Multiplicity TTV Period [days] TTV Amp [min] 


K00013.01 

1.76359 + 0.000001 

4633 

23 + 16.0 

1 



K00151.01 

13.4472 + 0.000047 

1426 

6 + 2.2 

1 



K00256.01 

1.37868 + 0.000002 

26406 

24 + 12.0* 

1 

47.5^* 

0.44^* 

K00258.01 

4.15745 + 0.000006 

1136 

5.1 +2.3 

1 



K00698.01 

12.7187 + 0.000006 

7514 

10.87 + 0.29 

1 



K00779.01 

10.406 + 0.000007 

14610 

12.1 + 5.1 

1 



K00823.01 

1.02841 +0.000004 

5867 

8.7 + 3.8 

1 

184'* 

0.85'* 

K00912.01 

10.8485 + 0.000031 

1651 

2.3 + 0.29 

2 



K01003.01 

8.36058 + 0.000005 

19878 

10.8 + 6.1 

1 



K01285.01 

0.937444 + 0.000003 

4971 

7 + 2.3 

1 

_AC 

_AC 

K01459.01 

0.692007 + 0.000001 

3876 

4.12 + 0.29 

1 

_AC 

_AC 

K01546.01 

0.917547 + 0.000002 

14150 

9.5 + 2.9 

1 

_AC 

_AC 

K01561.01 

9.08594 + 0.00004 

2184 

26 + 16.0 

1 



KOI 799.01 

1.73108 + 0.000002 

12542 

49 + 31.0 

1 



K02052.01 

4.22464 + 0.000028 

1026 

4.1 +2.0 

1 



K02867.02 

18.9407 + 0.0014 

14 

0.36 + 0.15 

2 




Table 1. Planet parameters. Values used from the Kepler Objects of Interest table from MAST. ' ^Mazeh et al.|(2oT3) . ^Szabo et al.|j2013) . ^ No value was 
given in the respective paper, but transit timing variations were still found to be significant. 


There is the obvious possibility that we introduced artificial 
signals in the light curves by cutting out transits. Any kind of 
Fourier analysis such as the Lomb-Scargle periodogram is suscep¬ 
tible to aliasing. To investigate this issue, we applied our transit cut¬ 
ting methods twice, once at the actual orbital frequency and once 
at a frequency that is an irrational multiple of the orbital frequency, 
Vinjected = 1.3575656 Vopbit- Thls value was chosen such that nei¬ 
ther Vinjectcd nor any of its harmonics fall near the orbital frequency 
or its harmonic, thus obscuring any artificial signals that we were 
attempting to create. The light curves were then prepared (normal¬ 
ized and cut) in the same fashion as for the original analysis. From 
this point, the periodograms were created with the same methods 
described previously, and were examined for newly created signals. 
In the power spectra, we were unable to detect artificial signals cre¬ 
ated by our cutting method near the frequency Vi„jectcd or any multi¬ 
ple thereof. Thus, we can rule out our data preparation methods for 
the origins of the higher order harmonic. 

We investigated one more possibility that could lead to arti¬ 
facts in the periodograms. A potential issue might arise for sys¬ 
tems of high eccentricity. In such systems, secondary transits can 
be shifted significantly from where we would expect them for a 
more circular orbit. This would cause the transit to bleed over and 
not be cut properly. With these extra data-points now contributing 
to the periodogram, higher order signals may be introduced. Sim¬ 
ilarly, transit timing variation could cause the same issue. We ran 
several tests with larger transit windows. We ranged our transit du¬ 
rations from a factor of 1.5 to a factor of 3 of the transit width. In 
all of these tests, the higher order harmonics did not disappear, but 
expectedly the signal to noise ratio was lower. 

Not all signatures in the light curves may necessarily be si¬ 


nusoidal. We tested generating synthetic light curve data contain¬ 
ing both sinusoidal and non-sinusoidal periodic signals. From this 
test, higher order harmonics can be recovered with varying results. 
As expected, the shape of the non-sinusoidal signal determines the 
relative powers of the higher order harmonics. Due to the nature 
of these signals, amplitude measurements for ellipsoidal variation, 
and beaming can be significantly alfected. These effects are impor¬ 
tant as they have recently been used to determine masses of planets. 
We note that in |Shporer et al.|p014[ >, beaming measurements were 
found not to agree with radial velocity measurements for mass es¬ 
timates. 


4 TIDAL EFFECTS 

If the third harmonic was entirely tidal effects, it would be expected 
to be similar to the effect on a single star in a binary system as 
described by |Morris & Naftil^ ( |1993| l: 

flsc = -25^*(2 -I- 1/32(3 -;/*) (1) 

Where is the cosine term of the harmonic on triple the frequency 
of the planet, /r* is the linear limb darkening coefficient, A is the 
semi-major axis of the planet, i is the inclination of the planet’s 
orbit, Mp and M* are the masses for the planet and star, t, is the 
gravity darkening coefficient, and is the radius of the star. Us¬ 
ing a least-squares fit with a three-harmonic model, we measured 
the amplitude of the third harmonic. In |Morris & NaftilanH1993] l 
tidal effects between binary stars were expected to produce high 
order harmonics at frequencies that are multiples of the orbital pe¬ 
riod. These tidal effects were applied to Kepler-13b in|Shporer et al.| 
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KOI Name 

Expected a^c 

Fitted a^c 

K00013.01'’' 

2.61 ■ 10~’ 

2.39 ■ 10“'’ 

K00151.01 

1.33 ■ 10~“ 

5.95 ■ 10“'’ 

K00256.01® 

- 

6.00 ■ 10“^ 

K00258.01 

1.03 ■ 10-^ 

8.90 ■ 10“^ 

K00698.01 

3.24- 10““ 

7.20 ■ 10“'’ 

K00779.01 

1.90- 10“^ 

2.65 ■ 10“'^ 

K00823.01 

8.03 ■ 10“^ 

9.88 ■ 10“'’ 

K00912.01 

2.55 ■ 10““ 

3.24 ■ 10“^ 

KOI 003.01 

3.27 ■ 10“^ 

5.10- 10“^ 

KOI 285.01 

1.16- 10“* 

2.06 ■ 10“^ 

K01459.01 

2.44 ■ 10“* 

1.47- 10“^ 

K01546.01 

4.88 ■ 10“’' 

3.00 ■ 10“'* 

K01561.01 

1.97 ■ 10“* 

7.52 ■ 10“^ 

K01799.01 

1.04- 10“^ 

5.95 ■ 10“^ 

K02052.01 

7.92 ■ 10“‘5 

8.10- 10“3 

K02867.02 

2.93 ■ 10“*'’ 

2.01 ■ 10“^ 


Table 2. The amplitudes for both the expected and fitted term of the cosine 
term for the third harmonic are in the form of relative flux. The expected 
amplitudes were calculated with the parameters from the Kepler Catalogue. 
The mass of many of these planets is unknown, so to approximate the mass 
ratios, the radii were used with the assuming that the star and planet would 
have similar density. All of the errors in the expected values are of the same 
order as the value. ^ The value quoted is for consistency with the other 
planets. If the mass ratio from |Shporer et alT]j2014) is used, a similar factor 
of 2 is found. ® The KOI-256 system contains a white dwarf, for which the 
limb darkening coefficients are unknown and assumed to be zero. 

( |2014^ , to describe the high order harmonics found. The third order 
signal was found to have twice the amplitude expected in the |Morri^ 
|& Naftil^jl993| l model. For the systems analyzed in this paper, all 
overshoot the expected amplitude. This brings up two possibilities: 
tidal effects between a planet-star and star-star are different enough 
to warrant further analysis, or there is an additional unknown as- 
trophysical effect. It is unlikely to be entirely tidal, as many of the 
signals are many orders of magnitude larger than expected, and un¬ 
correlated with the value expected from theory. 


5 CONCLUSIONS 

We analyzed the light curves of 16 Kepler systems. Through the 
use of simple Lomb-Scargle periodigrams, we were able to find ev¬ 
idence for the beaming effect, ellipsoidal variations and reflection, 
as previously extracted successfully by other groups. In addition, 
of the 16 chosen planets, 13 exhibited higher order harmonics in 
the light curve, on frequencies that are multiples of the orbital fre¬ 
quency for which there is currently no explanation. Through careful 
testing, we conclude that the higher order harmonics found do not 
originate from our own data processing. Thus, we are left with the 
following possible explanations: 


Tidal effects. By applying the analysis in [Morris & Naftilan| 
( |1993t , we conclude that tidal effects cannot be the sole contribu¬ 
tor to high order harmonics. The measured values of the high order 
harmonics exceeds what tidal effects predict by a range of a fac¬ 
tor of 2 up to many orders of magnitude. It is possible that tidal 
effects between planet-star and star-star are significantly different, 
but unlikely to make up for the huge orders of magnitude. 

Non-sinusoidal periodic light curve variations. In our tests, 
we found that non-sinusoidal variations can indeed replicate large 
amplitude higher harmonics, and depending on the shape of the 
variations, different combinations of higher order harmonics can be 
excited. Non-sinusoidal light curves might for example be caused 
by non-isotropic reflection or thermal emission from the planet’s 
surface. 

An error in the Kepler pipeline Although unlikely, there 
could in principle be an error introduced early on in the Kepler 
pipeline. A more sophisticated pixel-level analysis of Kepler data 
could provide further insight into the signal’s origin. 

At the present, we speculate that the most likely cause of these 
high order harmonics are non-sinusoidal periodic light curve vari¬ 
ations. We note that there is a weak correlation between the am¬ 
plitude of the high order harmonics and the planet’s radius, but no 
obvious correlation to period. If the physical origin of these har¬ 
monics is fully understood, we are likely to gain new insight into 
the planetary properties. Until then however, one should exercise 
caution when determining planetary properties from the shape of 
the light curve as the harmonics of unknown origin affect the mod¬ 
elling of known astrophysical phenomena. 
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